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We theoretically investigate laser induced quantum transport in a two-level quantum dot attached 
to electric contacts. Our approach, based on nonequilibrium Green function technique, allows to 
include thermal effects on the photon-induced quantum transport and excitonic coherent dynam- 
ics. By solving a set of coupled integrodifferential equations, involving correlation and propagator 
functions, we obtain the photocurrent and the dot occupations as a function of time. The character- 
istic coherent Rabi oscillations are found in both occupations and photocurrent, with two distinct 
sources of decoherence: incoherent tunneling and thermal fluctuations. In particular, for increasing 
temperature the dot becomes more thermally occupied which shrinks the amplitude of the Rabi os- 
cillations, due to Pauli blockade. Finally, due to the interplay between photon and thermal induced 
electron populations, the photocurrent can switch sign as time evolves and its stationary value can 
be maximized by tunning the laser intensity. 



I. INTRODUCTION 

Quantum transport in semiconductor quantum dots 
and molecular systems is a subject of intense 
study nowadays. 1 These nanoscalcd devices provide a 
formidable environment to study fundamental aspects of 
quantum physics, involving many-body correlations and 
light-matter interaction in regimes out of equilibrium. 2 
These systems have a great potential to form a new gener- 
ation of optoelectronic devices based on the unique elec- 
tronic structure that arises from the quantum confine- 
ment. For instance, quantum dots can produce a wealth 
of visible colors depending upon its size, even white light 
with relatively high efficiency 3 and potential to intcgra- 
bility with nanoelectronics. 4 Additionally, with the great 
technological advances in the manufacturing of semicon- 
ductor quantum dot system, it became possible to coher- 
ent monitor and control electron populations in two-level 
systems via different pump-probe techniques. 5-8 In all 
these experiments the main signature of quantum coher- 
ent nonlincarity is Rabi oscillations, which has no clas- 
sical counterpart. More recently, Rabi oscillations was 
also reported in organic light-emitting diode. 9 Such co- 
herent optical manipulations constitute a fundamental 
ingredient to quantum information processing in solid 
state devices that use electron-spin or excitonic states as 
qubits. 10 Interestingly, holes in semiconductor quantum 
dots have been revealed as an alternative to electrons in 
the manufacturing of spin qubits. 1 

It was originally demonstrated by Zrenner et ai. 12 that 
coherent Rabi oscillations in a two-level quantum dot 
photodiode can be monitored by photocurrcnts. Addi- 
tionally, it was proposed that a photocurrent in a self- 
assembled quantum dot photodiode can become spin- 
polarized due to an effective exchange interaction via 
biexciton state. 13 This result points out the potential- 
ity of the present system to future spintronic devices. It 
was also observed that the double dot structures present 
the ability to increase the coherence time of indirect 
excitons. 14 Recently, thermal effects on the excitonic 
Rabi rotations in a quantum dot system were investigated 



experimentally. 15 ' 16 It was evidenced acoustic phonons as 
the main source of damping of the Rabi oscillations. 

In the present work we analyze how the temperature 
of nearby contacts tunnel coupled to a two-level quan- 
tum dot affects the coherent optical dynamics. Apply- 
ing nonequilibrium Green function technique to a mi- 
croscopic Hamiltonian model, we write a set of coupled 
integrodifferential equations that describes the coherent 
evolution of the electron- hole populations in the dot. A 
resonant laser field drives the electron-hole dynamics and 
generates a photocurrent. We find two contributions to 
the current. The first one comes from electrons in the 
dot that tunnel to a contact. This current is positive and 
is mainly induced by the laser field. The second current 
component is related to electrons in the reservoir that ac- 
quires enough thermal energy to tunnel into the dot. This 
second contribution charges the dot, thus generating ad- 
ditional features not yet reported in the literature. Our 
main findings include the suppression of the Rabi oscilla- 
tion as the temperature of the nearby contact increases, 
a negative photocurrent due to a backwards charge flow, 
and a maximum photocurrent value achieved when the 
laser intensity is comparable to the mismatch between 
the excited dot level and the contact chemical potential. 



II. MODEL AND FORMULATION 

Figure (1) illustrates the system considered. It is com- 
posed of a quantum dot attached to a left and to a right 
electron reservoirs in the presence of a source-drain bias 
voltage. A laser field shines the dot, thus generating 
electron-hole pairs in it. The electrons in the conduction 
band and the holes in the valence band can tunnel out 
the quantum dot to the left and to the right electrodes, 
respectively. This results in a photocurrent throughout 
the system. 

The modeling Hamiltonian is given (per spin) by H a = 
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FIG. 1. (Color online) Sketch of the system studied. A quan- 
tum dot is tunnel coupled to both left and right reservoirs. 
The left reservoir illustrated has chemical potential and 
temperature T. Due to the proximity of the conduction band 
level £2 to fi2, the dot population thermally fluctuates. This 
induces nonlinearities in the photocurrent driven by a laser 
field with resonant energy hw. 

H Dt(T + H T ,a + iih.a + H Rt<r , with 

where ti a is the dot level for spin a in the valence (i = 1) 
or the conduction (i = 2) band. The operators di„ (d\ a ) 
annihilates (creates) one electron in level i with spin a. 
The parameter 7 gives the optical transition between va- 
lence band and conduction band in the quantum dot. 
This parameter can be controlled by the intensity of the 
incident radiation. To couple the dot to fcrmionic reser- 
voirs we use the tunneling Hamiltonian 

E ?,° = E Z> c L** + f4^). ( 2 ) 

% ki 

where Ck il7 (ct.-) annihilates (creates) one electron in the 
right (i = 1) or the left (i = 2) lead. 17 The parameter ti 
gives the dot-leads coupling strength. Finally, the free- 
electron energies of the electrons in both leads are given 
by 
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H L ,<J + Hr !<7 = ^ 2J e fciCr C l iC r C fciCr- ( 3 ) 

»=1 ki 

In the present model we assume that the tunneling 
rates are larger than spontaneous emission rates, so that 
electron-hole recombination will be neglected. 

Our main task is to explore the effects of the reser- 
voirs temperature on the Rabi oscillations and the pho- 
tocurrent. To this goal we must find the lesser and 
retarded Green functions of the quantum dot, i.e., 

G< CT (M) = i(4„(f)(k(*)> an d G^ftf) = -iO(t - 
t')({dicr(t),djcr(t')}). Note that the occupations of the 
levels e, CT are given by rii a (t) — \m{Gf icr (t, t)}, while 
the photocurrent is defined as 7j CT = —e(Ni a ) = 



-eidH^^ig]), (h = 1) with N ia = 4 iCT c fei(T being 
the total number of particles operator. Following Rcf. [2] 
one can show that 

I 2a (t) = -eT 2(T n ia {t) + 2eRe{$ r 22tT } (4) 
where <5> r 22ty {t) = J^dhG^^t^fo^t^t) and 

faitut) = iT ia J ^/ i ( e )e- <e ( tl - t ), (5) 

with i = 1,2. Here /j(e) is the Fermi function to i- 
th reservoir and = 2ir\ti\ 2 pi a is the tunneling rate 
with pi a being the density of states of the corresponding 
reservoir for spin component a. The present formalism 
allows the inclusion of ferromagnetic leads by considering 
spin-dependent tunneling rates r^. 23 According to Eq. 
(4) the current at time t has two contributions, one that 
is instantaneous and proportional to the dot occupation 
ni a (t) (I ut) and a second one that involves the whole 
history of the system (ij n ). In this second term, a time 
integral of the correlation function G 22a .(t,ti) weighted 
by a thermal dependent function 4>2o-(ti,t) should be 
carried on, ranging from — 00 to the present time. All 
the thermal effects arise via this memory integral. This 
contrasts to the density matrix approach used in quan- 
tum optics formulation that in general does not account 
for thermal and memory effects in the standard Markov 
approximation. 24 ' 25 

Calculating the time derivative of Gfj a (t,t) we arrive 

at 

i^G<(t,t) = M ff (i)G<(M) - <Mi), (6) 

where the lesser Green function is written in a vector- 
like form G< = [G< G< a , G< „, G< CT ] T and * a = 

[*ii<r. *'i2„> *^F- Hcre = *«(*) ~ 

with 

$T.(t)= f dM%,,(Mi)M*i.*)> (7) 

J —00 

and 

$?■(*)= f dh^tJGf^t). (8) 

J —OO 

The matrix in Eq. (6) is given by 



M„(t) = 





— rye~ iLJt 
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with uiij — e, — ej and r CT = + r2cr. It is yet valid 
to note that in the absence of the reservoirs, Eq. (6) re- 
duces to the well known semiconductor Bloch equations. 2 
In order to determine $<r(i) we need the retarded and 
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advanced Green functions G[j a (t,t'). Taking the time 
derivative with respect to t' we obtain 



■P a (t')G r a (t,t'), (10) 



where G r a = [G r lla , G r 12a , G r 21a , G r 22a \ T , x+ and x~ are 
the two-component Pauli spinors, x+ = [1; 0] T and X- = 
[0, 1] T , and the matrix ~P a (t') is defined according to 



P«r(t) = 



7 e 



5i 7e 
iut' g. 



iojt' 



7e 






Ji 7e 



iujt 
2 



(11) 



with (5; = e; — ^r^. Solving numerically Eqs. (6) and (10) 
we obtain the occupation n^t) and the photocurrcnt. In 
what follows we present our results. 



III. PARAMETERS 

In order to keep the generality of our results, we ex- 
press the time in units of to = H/Tq, where To is the 
tunneling rate between dot and reservoirs. For simplic- 
ity we assume the wideband limit, where the tunneling 
rates are energy independent and we set = T 2 a = IV 
The current unit is given by Iq = eYo/h and all en- 
ergies are in units of IV Experimentally, Tq depends 
on the tunnel barrier and it can be easily controlled 
by an external gate voltage. We find for quantum dot 
systems To ~ lfieV — 100/ieV, 18-21 which results in 
Io ~ 0.24nA — 24nA. The time to ranges in the interval 
to ~ 6.5ps (T = lOOfieV) - 0.65ns (r = lfxeV). Addi- 
tionally, the quantum dot valence and conduction band 
levels are given by e\ c = e\ = — lOOFo and £2,7 = £2 = 
2IV respectively. 22 Both levels are measure with respect 
to the chemical potential fj, 2 = 0, which is taken as our 
energy reference. Finally, we assume fc^T ~ O.lTo — 2Fo 
and 7 = 7Fo, except when those parameters explicitly 
change in the plots. 



IV. RESULTS 

Figure 2(a)-(b) shows the evolution of the electron and 
hole occupations in the quantum dot for differing tem- 
peratures UbT. At t = the valence band level is fully 
occupied with n\ = 1 while the conduction level is little 
populated with n 2 w 0.1. This small occupation comes 
from the proximity of the level e 2 to the Fermi energy of 
the left reservoir, which allows thermal excited electrons 
to tunnel into the dot. Initially (t — 0) the quantum dot 
occupation is calculated through 



n-«r(t = 0) 



J^-OUe), (12) 



where the lesser Green function is given by the Keldysh 
equation G< a (e) = G[ iCT (e)S< (e)G^(e), where G^ is 
the retarded (advanced) Green function of the dot at- 
tached to the leads without laser field and Y,f a = iT ia fi. 
As UbT increases, electrons in the left electrode acquire 
enough thermal energy to enhance the population n 2 at 
t = 0, while ni remains the same due to ej <C £f- When 
the system starts to evolve in the presence of a laser field, 
the occupations ni and n 2 develop the characteristic Rabi 
oscillations. In the small temperature regime these oscil- 
lations are more pronounced for small times and become 
suppressed as the time increases. This is due to the deco- 
herence imposed by the tunneling between dot and reser- 
voirs. 

As the temperature increases, the amplitude of the 
Rabi oscillations shrinks for all times. This is related to 
the enhancement of n 2 with temperature. With the level 
e 2 becoming more populated, the Pauli exclusion princi- 
ple makes it more difficult to one electron with the same 
spin in the valence band (ei) to jump to the conduction 
band (£2)- So we observe two sources of suppression to 
the coherent Rabi oscillations: (i) coupling to reservoirs 
and (ii) thermal activated Pauli blockade. 

The photocurrent seen in Fig. 2(c), at least to some 
extent, reflects the n 2 behavior. It oscillates in time with 
a decreasing amplitude, tending to a stationary nonzero 
value. Additionally, the amplitude of the Rabi oscil- 
lations is also reduced as fc^T increases, following the 
behavior of n 2 . Interestingly, for small enough temper- 
atures and shorter times, the photocurrent oscillations 
attain negative values, which corresponds to an unusual 
flow of electrons from the reservoir into the dot (see solid 
line, ksT = 0.1 To, around t = 0.5to). In order to gain 
further insight of this effect, in the inset of Fig. 2(c) 
we show separately the current components I ou t and Jj n . 
The outgoing current is positive which means that elec- 
trons are flowing from the dot to the reservoir. The in- 
coming current gives a negative contribution to the cur- 
rent, which corresponds to electrons flowing in the op- 
posite way, i.e., from the electrode into the dot. Around 
t = 0.5to, the Ii n component presents a dip, which pulls 
down the total photocurrent, making it negative. When 
ksT increases, this dip is suppressed and the photocur- 
rent assumes positive values for all times. 

Figure 3(a) shows how n\ and n 2 evolves func- 
tion of the parameter 7, for differing temperatures. For 
all fcgT values we observe n\ = 1 for 7 = 0, while n 2 in- 
creases with ksT for 7 = 0. This enhancement of n 2 with 
temperature comes from the thermal excited electrons in 
the reservoir that acquires enough energy to jump into 
the dot as kgT increases. Both n\ and n 2 are obtained 
via Eq. (12) for 7 = 0. In the presence of the laser field, 
n 2 increases monotonically with 7, while n\ is initially 
suppressed and then it is enhanced, thus developing a 
minimum around 7 f=a e 2 — [i 2 = 2To- Note also that n 2 
presents a further enhancement near 7 = 2Fo. In the in- 
set of Fig. 3(a) we show the sum n — n\ + n 2 against 7. 
When the resonant condition 7 = e 2 — fi 2 is attained, the 
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FIG. 2. (Color online) (a)- (b) Occupation of the levels ei and 
£2 and (c) photocurrent as a function of time for differing 
temperatures. As time evolves all these quantities exhibit co- 
herent Rabi oscillations. The amplitude decreases with time 
due to the incoherent tunneling between dot and electrodes. 
The oscillation amplitude is also suppressed as ksT increases 
due to a thermal induced Pauli blockade. For t « 0.5to the 
photocurrent is dominated by a backwards current, thus be- 
coming negative. In the inset we show separately the in and 
out current components. 



total occupation presents a steeper enhancement for low 
temperature. For larger fc^T a more broaden increasing 
is found. The occupation profiles will be more clearly un- 
derstood looking at the current components in the next 
plot. 

In Fig. 3(b) we plot separately the current components 
lout, lin and the total current I — I ou t + hn- While 
the outgoing current follows n^, the incoming current 
is strongly increased (in modulus) around 7 ps £2 — M2- 
As a result, the photocurrent is suppressed due to this 
backwards current, thus developing a peak close to 7 = 
(2 — A*2- Increasing even further the temperature, the 
thermal fluctuations of the reservoirs yield to a more ef- 
fective injection of electrons into the dot. This makes 
Ii n starts at higher absolute values for 7 < 2To. This 
amplification of Jj n suppresses the photocurrent when 
compared to its low temperature profile. In the presence 
of a laser field in resonance with the difference £2 — ei, 
doublets emerge in the spectrum of the system, 21,26,27 as 
illustrated in the drawn of Fig. 3(b). 28 As the laser in- 
tensity increases, the lower energy peak of the doublet 
eventually attain resonance with the reservoir chemical 
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FIG. 3. (Color online) (a) Occupations ni and ri2 and (b) 
current components against 7. While 712 increases monoton- 
ically with 7, m is initially suppressed, reaches a minimum 
and then increases slightly. Oppositely, the photocurrent in- 
creases with 7, reaches a maximum and then becomes slightly 
suppressed. This nonmonotonic behavior in both n\ and I can 
be understood looking at the current component I in . When 
7 ~ £2 — /J, = 2ro, electrons in the reservoir can tunnel into 
the dot, thus populating £2. This enhancement (in modu- 
lus) of the backwards current suppresses the photocurrent for 
7 > 2Fo and increases further ri2. In the inset we plot ru+na- 
Note that the dot charges for 7 > 2Fo. 

potential at 7 = £2 — [ii- This allows electrons to reso- 
nantly tunnel from the lead into the dot, thus generating 
a backwards current that suppresses the total photocur- 
rent and increases the 712 population. When the lower 
peak lies below /i2, the enhancement of fcgT tends to de- 
populate this channel, consequently suppressing Zj„, as 
seen in Fig. 3(b) for 7 > 2r . 

Finally, it is valid to point out that the Ii n current com- 
ponent plays a role in the transport whenever £2—7 < M2j 
which allows electrons in the reservoir to resonant tun- 
nel to the dot. It is possible to entirely suppress the 
incoming current by moving £2 high enough above /i2, 
so that £2 — 7 > M2- For this regime electrons can flow 
only in one direction, i.e., from the dot to the reservoir, 
thus reducing the backwards charge flow. Fig. 4 shows 
the photocurrent against time for different £2 values. For 
£2 = 2Fo and 5Fo the channel £2 — 7 lies below /J2, as 
illustrated in the energy diagram at the lower part of the 
panel. This results in a relatively high component 
(see the inset). On the contrary, for £2 = lOTo and 20Fo 
we find £2—7 higher than /J2 (see the upper energy di- 
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agram sketched), which suppresses |7j„| and makes the 
total current larger. 
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FIG. 4. (Color online) Photocurrent against time for different 
£2 values, with respect to ^2. As £2 increases the photocurrent 
is amplified. This is due to the suppression of the incoming 
current for larger £2. The temperature adopted is fcsT = 
O.lTo. In the inset we show the incoming current component 
for all the £2 values used. As £2 enlarges, the incoming current 
tends to zero. 



V. CONCLUSION 



In conclusion, via nonequilibrium Green function tech- 
nique we have investigated the coherent dynamics of 
electron-hole pairs in a quantum dot tunnel coupled to 
Fcrmionic reservoirs. As temperature increases, the ther- 
mal excited carriers in the left reservoir acquires enough 
energy to tunnel into the dot. This gives rise to an en- 
hancement of electronic dot population, which results in a 
thermal activated Pauli blockade that suppresses slightly 
the Rabi oscillations. This effect is strongly dependent 
on the temperature of the reservoirs and on the mismatch 
between e 2 and /Z2. Finally, a nonlincarity signature is 
found in the current against 7. This results from a dou- 
blet that brings into resonance a transport channel with 
the chemical potential fi 2 - This laser induced resonance 
generates a competition between outgoing and incoming 
currents in the quantum dot that yields to the observed 
nonlinearities. 

The authors acknowledge the Brazilian agencies 
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